Effects of interaction on an adiabatic quantum electron pump 
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We study the effects of inter-electron interactions on the charge pumped through an adiabatic 
quantum electron pump. The pumping is through a system of barriers, whose heights are deformed 
adiabatically. (Weak) interaction effects are introduced through a renormalisation group flow of 
the scattering matrices and the pumped charge is shown to always approach a quantised value at 
low temperatures or long length scales. The maximum value of the pumped charge is set by the 
number of barriers and is given by Q ma x = — 1. The correlation between the transmission and 
the charge pumped is studied by seeing how much of the transmission is enclosed by the pumping 
contour. The (integer) value of the pumped charge at low temperatures is determined by the number 
of transmission maxima enclosed by the pumping contour. The dissipation at finite temperatures 
leading to the non-quantised values of the pumped charge scales as a power law with the temperature 
(Q — Qint oc T 2a ), or with the system size (Q — Q ln t c< Lj 2a ), where a is a measure of the interactions 
and vanishes at T = (L s = oo). For a double barrier system, our result agrees with the quantisation 
of pumped charge seen in Luttinger liquids. 



PACS number: 73.23.Hk, 72.10.Bg, 73.40.Ei, 71.10.Pm 
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In the last few years, there has been a lot of interest 
in the phenomenon of electron transport at zero bias via 
the mechanism of adiabatic pumping. The parameters 
of a system are slowly varied as a function of time, such 
that they return to themselves after each cycle, and the 
net result is charge transport [1] . The idea of quantised 
charge pumping was first introduced by Thoulcss [2] for 
a spatially periodic system, where the quantisation could 
be proved using topological arguments. This was further 
investigated theoretically in Refs. [3]. 

Experimentally, the possibility of transferring an in- 
teger number of charges through an unbiased system 
has been achieved in small semi-conductor dots in the 
Coulomb blockade (CB) regime [4,5]. It has also been 
achieved in surface-acoustic- wave based devices [6], but 
here too, the quantisation has been attributed to CB. An 
adiabatic electron pump that produces a dc voltage, in 
response to a cyclic change in the confining potential of 
an open quantum dot has also been experimentally con- 
structed [7]; however, here the charge is not quantised. 

The open dot quantum pumps, where the quantum in- 
terference of the wave- function, rather than the Coulomb 
blockade, plays a dominant role, has been of recent the- 
oretical interest. A scattering approach to such a pump 
was pioneered by Brouwer [8] and others [9,10], who, 
building on earlier results [11], related the pumped cur- 
rent to derivatives of the the scattering matrix. Using 
this approach, there have been several recent theoretical 
studies [12-31]. The conditions under which the charge 
is (almost) quantised in open quantum dots, has been ex- 
plored in detail in Refs. [32-35]. Ref. [32] shows that the 
pumped charge is produced by both nondissipative and 
dissipative currents. It also shows that charge quantisa- 
tion is achieved when the dissipative conductance van- 
ishes. This is also emphasized in Ref. [35], where the 



adiabatic charge transport in a Luttinger liquid has been 
studied. In Refs [33,34], the correlation between the con- 
ditions for resonant transmission and quantised charge 
has been emphasized. 

In this paper, we include the effects of inter-electron 
interactions along with the interference effects, using a 
renormalisation group (RG) method developed recently 
[36-39] for the flow of the scattering matrices. We show 
that at low temperatures, the inclusion of interactions, 
leads to quantisation of the pumped charge. As seen ear- 
lier in Refs. [32] and [35], we find that the charge pumped 
is a sum of two terms - a quantised term, and a dissipative 
term. The quantised term, which is topological in nature 
and is due to quantum interference is always present. 
For non- interacting electrons, the dissipative term is also 
large and masks the quantisation due to the first term. 
However, the RG flow of the S-matrices, due to the in- 
teraction, causes the dissipative term to vanish at long 
length scales or low temperatures making the quantisa- 
tion of the pumped charge visible. 




FIG. 1. Schematic diagram of a multiple dot system (rid 
dots or rib + 1 barriers) defined on a two dimensional electron 
gas. The barriers forming the dots are denoted as Vt and the 
gate voltages controlling the density in the dots are denoted 
as Vq. 

The quantum pump is formed by a system of coupled 
quantum dots (shown in Fig.l), where the barriers form- 
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ing the dot are periodically modulated as 

Vi = V\ = Vo + V p cos(wt), i < rib/ 2 for rib = even, 

i < rib/2 + 1 for rib = odd, 
Vi = V2 = Vo + V p cos(uit + <5), for the remaining. (1) 

Clearly, for rib barriers, we have rid = rib — 1 dots. Here 
lo is related to the time period as u = t/2tt and S is the 
phase difference between the two time- varying potentials. 
Since this potential breaks the parity symmetry, it allows 
the shape of all the dots to be varied. As has been empha- 
sized earlier, this is a necessary condition to get non-zero 
pumped charge. 

Although we use S- function barriers, we expect the cal- 
culations to be robust to changing the form of the barri- 
ers. Following the work of Ref. [40], we also expect these 
results to be robust to weak disorder. We treat the dot 
within a one-dimensional effective Hamiltonian [41], since 
they are coupled to the leads by a single channel quantum 
point contact. (Although in real quantum dots, transmis- 
sion properties can fluctuate and energy level spacings 
can vary and can even be chaotic, this modelling works 
to understand simple physical features of higher dimen- 
sional quantum dots.) The width of the dots (effectively 
the width of the quantum well that we use) is given by a. 
We are mainly interested in the region where Vb < Ep 
since we wish to study the resonant tunneling regime and 
not the CB regime. 

The effective two-channel S-matrix for this system of 
rib barriers can be written as 

s _(re» fe**\ 

where the parameters r, /, 9, 9' and (j) are functions of the 
Fermi energy Ep and the amplitudes of the time- vary- 
ing potentials Vi{t). Their explicit forms can be found, 
in terms of the parameters of a single well, (in the adia- 
batic limit), by solving the time-independent Schrodingcr 
equation for the potential Vi(t) given in Eq. 1, for each 
value of t. The reflection coefficients are not the same 
(phase can differ) because the time- varying potentials ex- 
plicitly violate parity. The potential also violates time- 
reversal invariance. But since in the adiabatic approxi- 
mation, we are only interested in snapshots, at each value 
of the time, the Hamiltonian is time-reversal invariant 
and hence, the transmission amplitudes are the same for 
the '12' and '21' elements in the S-matrix. 

By the Brouwer formula [8], the charge pumped can 
directly be computed from the parametric derivatives of 
the S-matrix. For a single channel, it is given by 



where Sij denote the matrix elements of the S'-matrix and 
Vi and V2 are the time derivatives of the V\ , V2 given in 
Eq. 1. Thus, the pumped charge is directly related to 
the amplitudes and phases that appear in the scattering 
matrix. For the form of the S'-matrix given in Eq. 2, this 
is just 

Q = ^j\o-f 2 (o~4>)\dt (4) 

where the first term is clearly quantised since e has to 
return to itself at the end of the period. So the only 
possible change in 9 in a period can be in integral multi- 
ples of 2n -i.e., 9(t) — > 9(0) + 2nn. The second term is 
the 'dissipative' term which prevents the perfect quanti- 
sation. It is easy to see the analogy of Eq. 4 with Eq. 19 
of Ref. [32]. 

The correlation between resonant transmission and the 
pumped charge has been studied earlier by several groups 
[16,26,30,33,34]. A pumping contour [33,34] can be de- 
fined in the plane of the parameters, V\ and V2, as the 
closed curve traced out in one cycle of modulation of the 
barriers (t — > t + r). As we shall see below explicitly, 
the value of the first term in Eq. 4 is set by the number 
of transmission maxima enclosed by the pumping con- 
tour, whereas the second term depends on how peaked 
the transmission function is, around its maxima. Wc 
find that the more peaked the transmission function, the 
smaller is the dissipation. 

The aim of this paper is to study the effect of electron- 
electron interactions on the charge pumped. In Ref. [35], 
a bosonised approach was used; here, the barriers have 
to be either in the weak back-scattering or the strong 
back-scattering limit. But for time-dependent potentials, 
within a single pumping cycle, the barriers can range 
from the weak barrier to strong barrier limits, if the 
perturbing potential is large. Contours with such large 
perturbing potentials are required for enclosing multi- 
ple transmission maxima, which, in turn, lead to larger 
pumped charges. Hence, it is much more convenient to 
introduce interactions (provided they are weak) pertur- 
batively and directly obtain the renormalisation group 
(RG) equations for the entries of the S-matrix. This 
method has been used extensively in Refs. [37,39] and we 
directly present the RG equations for the effective two- 
channel S-matrix given in Eq. 2 as follows - 

^ = M - SM^S . (5) 
Here M is a diagonal 2x2 matrix given by 

M = ( 5« ir ? ) (6) 

and I is the length scale. The dimcnsionlcss constants 
ctj are small and positive and are related to the Fourier 
transform of the short-range density-density interaction 
potential Vi(x) in the two channels as 
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(Fi(o)- v;(2fcf) 
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(7) 



where we assume that the Fermi velocity is vp = Tifep/m. 
Thus, a.i is a measure of the inter-electron interactions. 
It can also be related to the bosonisation parameter K 
as 
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(8) 



When the barriers are symmetric (i.e., at some point in 
the pumping cycle when V\ = Vi ) ) , parity is not violated 
and we can set 9 = 9' in Eq. 2. We will always run 
the RG equations at this point. For the computation of 
the pumped charge, the origin in t is not relevant. We 
only need to compute the charge pumped in one cycle. 
It is also natural to expect the two channels to have the 
same inter-electron interactions; hence on = a. In that 
case, we obtain the following explicit equations for the 
RG flow of the transmission (and reflection) amplitudes 
and phases. We find 



dl 



(9) 
(10) 



for the flow of the transmission amplitude and phase and 



J = |r(l-r 2 -/ 2 cos2(^-0)) 



dO a „ 9 . „. , 



(ii) 

(12) 



for the flow of the reflection amplitude and phase. (Note 
that when 9^9', two similar sets of equations can be 
obtained, both for 9 and for 9' using S12 and Su for 
one and S22 and S21 for the other.) However, the condi- 
tion of unitarity on the S-matrix in Eq. 2 implies that 
4> — 9 = ir/2 + 2nir. This simplifies the RG flow equations 
for the reflection amplitude and phase, so that they be- 
come similar to those for the transmission amplitude and 
phase, i.e., 



dr 

— = ar(l 
dl 

d l = o 

dl 



(1-r 2 ), 



(13) 
(14) 



The entries of S therefore become functions of the length 
scale I or equivalcntly of L, where I = \n(L/d), and d is 
the short-distance cutoff which we take to be the average 
interparticle spacing. The RG flow can also be considered 
as a flow in the temperature since the length scale L can 
be converted to a temperature T using 
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as the thermal length. The interpretation of Lt is that 
it is the thermal length beyond which the phase of an 



electron wave-packet is uncorrelated. Thus Lt is the 
thermal phase coherence length. Hence, the RG flow 
has to be cutoff by either L T or the system size L s , 
whichever is smaller [39]. In all our numerical compu- 
tations, we start the RG flow at the microscopic length 
scale d, which is taken to be the inter-particle separa- 
tion. We then renormalise to larger and larger length 
scales, stopping at either Lt or L s , whichever is smaller. 
Note that for a length L = 100/zm, the equivalent tem- 
perature as computed from Eq. 15 is T = 6.5mK (for 
vp = 3 x 10 cm/ sec, which is the typical value of the 
Fermi velocity in GaAs). Both of these are experimen- 
tally achievable. Hence, whether L s or Lt is reached 
first depends on the given system. If Lt > L s , we can 
continue the RG flow until we reach L T , where phase 
coherence is lost due to thermal fluctuations. 

In the rest of the paper, we compute the transmissions, 
the phases and the quantised charges as a function of the 
various parameters in the theory. We then show how the 
elements of the S'-matrix, and consequently, the trans- 
missions, the phases and the pumped charge change as a 
function of the RG flow. Since the RG flow is also a flow 
in temperature, we show how the pumped charge changes 
as we decrease the temperature and reaches quantised 
values at very low temperatures. However, note that 
the RG flow is cutoff by the physical size of the sys- 
tem at a temperature Tl s — hvp/kBL s and one cannot 
go to lower temperatures. However, one can still con- 
tinue the RG flow upto the length scale Lt, until the 
system loses thermal phase coherence. ( For illustrative 
purposes, though, we sometimes take the limit that the 
system size goes to infinity. In that case, the tempera- 
ture can be lowered to T — 0.) We will see that as we 
change the temperature, the pumped charge shows power 
law behaviour as a function of the temperature until we 
reach the system size. Beyond that, the pumped charge 
shows power law behaviour as a function of the system 
size. In the next few subsections, we explicitly study the 
double barrier system n& = 2, the triple barrier system 
rib = 3 and the quadruple barrier system n& — 4. We 
then extrapolate the results to an arbitrary number of 
barriers. 

Strictly speaking, to remain within adiabatic approxi- 
mation under which the Brouwer formula is derived, the 
energy level spacing in the dots A has to be larger than 
the energy scale defined by the frequency of the time- 
varying parameter = Klu. It is only under this ap- 
proximation that the snapshot picture of studying the 
static S'-matrix for different time points within the pe- 
riod is valid. A better approach to go beyond the adia- 
batic approximation [24,42] is to use the Floquet states. 
Here we remain within the adiabatic approximation even 
though Ep > Vq, by choosing the spacing between the 
barriers a to be sufficiently small, so that the separation 
of the resonance level peaks (A > huo). 

• Single dot case or rib = 2 : 
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FIG. 2. Transmission F = f 2 as a function of the 
barrier strengths V\ and V2 for three different values of 
L = 1,100,1000 (Tl = 18if, 0.18K, ISmK) for the double 
barrier system. For Vo — 4.0, V p — 1.0, the pumping contour 
is shown as a dashed line. The other parameters are given 
by E F = 8.4, a = 0.3, uj = 1, o = 4 and 5 = ty/2. The 
pumped charges are 0.027, 0.29, 0.62 for the three values of L 
respectively. 



S = 



Ae 2la ^/B 
1/B 



1/B 

Ae -2iaVW/B 



(18) 



where A = ~i(V 1 + V 2 e 2iaVW ^)/2y/E^ -V 1 V 2 (e 2la ^ : - 
1)/4E F and B = 1 + i(V x + V 2 )/2y/E^ + V 1 V 2 {e 2ia ^ - 
l)/AEp and V\ and V 2 are the strengths of the two bar- 
riers. 
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FIG. 3. Q (in units of e) as a function of the RG length 
scale, L for three different values of the interaction strength 
a = 0.2, 0.25, 0.3 for the double barrier system. The values of 
the other relevant parameters are the same as in Fig. 2. The 
inset shows the change in the reflection phase in one pumping 
cycle to be 2n. 



Here, we compute the scattering matrix for two 6- 
function barriers at a distance a apart. To facilitate the 
generalisation to a larger number of barriers, we first ob- 
tain the M-matrix linking the incident wave to the left of 
a (5-function barrier at x — a to the transmitted wave to 
the right of the barrier, in terms of the strength V of the 
potential, the distance a and the energy of the incident 
wave Ep as [30] 
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iV_ „2iay/E^ 
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(16) 
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The S'-matrix which relates the outgoing waves to the 
incoming waves is simply obtained by rewriting the above 
matrix elements as 



S = 



M 21 /M u 
1/M n 



1/M n 
M 12 /M n 



(17) 



Since for a series of potential barriers, the transmit- 
ted wave can be simply obtained by multiplying the 
Af-matrices, it is now easy to obtain the S'-matrix for 
any number of barriers and thus identify the parameters 
r, /, 9, 9' and <fi in terms of the parameters of the potential 
scattering such as the distance a, the potential strengths 
Vi and the Fermi energy Ef- For two barriers, we find 
the elements of the S'-matrix as 
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FIG. 4. Scaling of the dissipative term Q — Qi nt as a func- 
tion of the length scale for the double barrier system. As seen 
in Eq. 20, for long length scales, the term that appears in 
the denominator can be ignored and the dissipation terms for 
different values of Ef fall on the scaling curve L~ 2a shown in 
the figure as a solid line. The contour chosen has Vo = 4.5, 
V p = 4.0 and S = tt/4. 

Following Ref. [16], to obtain numerical results, we set 
the width of the well a = 4 (in units of the inter-particle 
separation aj — 100 A") and u> = 1. We find, as expected 
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in the adiabatic approximation, that our results are in- 
dependent of uu and hence u> can be made as small as we 
wish. We note that our results are periodic in a for a 
fixed Ep (see Eq.16) and hence the width of the well can 
be increased. However, to remain within the adiabatic 
approximation, a is restricted by A > Klo. Also, for any 
value of a, one can find the maximum for the pumped 
charge by tuning Ep [16,30]. Our energy units are set 
by h = 2m = fcs = 1, where fee is the Boltzmann con- 
stant. So for GaAs, with a typical inter-particle spacing 
aj = 100yl o , and a typical effective mass m — 0.07m e , if 
we wish to set ai = 1, the energy unit is E = 5.6meV, 
(using E = h 2 /2maj) which corresponds to a temper- 
ature of T = 65°K (using E = k B T). The pumped 
charge for the double barrier system was obtained earlier 
in Refs. [16,30] and shown to be very small. Here, we 
shall see that the reason for this is that the transmission 
is not peaked about the transmission maxima; in fact, it 
is rather flat. Hence, there is a large dissipative term in 
the pumped charge. We shall also see how interactions 
make the transmission more peaked about the maxima, 
reduces the dissipation and enhances the possibility of 
quantisation. We compute the transmission as a func- 
tion of the barrier heights V\ and Vi and plot it in Fig. 
2, for a range of Vi between 2.0 and 6.0. The pumping 
contour [14] is defined as the closed path formed in the 
parameter space by the variation of the parameters in 
one cycle. Here, the varying parameters are V\ and V2. 
For this graphs, we choose the mean barrier height to be 
Vb = 4.0 and the perturbation to be V p = 1.0. We also 
choose 8 = 7r/2 to maximise the pumped charge. (For 
the triple barrier and quadruple barrier cases, different 
pumping contours are chosen depending on what we need 
to illustrate). These parameters, using Eq. 1, define the 
pumping contour, shown in the V\ — V2 plane in Fig. 2 
as a dashed line. The top surface in Fig. 2 denotes the 
transmission (at the Fermi energy Ep = 8.4, which is 
much above the barrier heights 4.0 ± 1.0), without any 
renormalisation of the barriers due to interactions. It is 
essentially flat and only a small part of the region of high 
transmission is enclosed by the pumping contour. Thus, 
although the pumping contour encloses the transmission 
maximum, which essentially implies that the first term 
in Eq. 4 is unity, (see inset in Fig. 3, which shows that 
the reflection phase change in one pumping cycle is 2n), 
the second dissipative term is very large and the pumped 
charge is vanishingly small (Q = 0.027). 

We now use the RG equations in Eqs.9-14. Since the 
phases do not change under the RG flow, we only need 
to use Eq. 9, to compute the change in the transmis- 
sion as a function of the length scale. We start the RG 
flow at L = 1 = d = 4 x 100A°, (a = 4), which is the 
width of the well, and then renormalise to longer length 
scales, until we reach the size of the sample, which is 
at best 100/iTO. At this length scale, the entire sample 
is phase coherent. In terms of temperature, the renor- 



malisation is from T = 18K to T = 6.5mK, which is 
well within the range of experimental feasibility. ( How- 
ever, in many graphs, we show a much larger range of the 
renormalisation, assuming that the samples can be made 
phase coherent, over much longer length scales, simply 
for illustrative purposes.) We find that as the RG flow 
proceeds, the transmission maxima get more and more 
peaked as seen in Fig. 2, where we have plotted the 
transmissions at different points along the RG flow. The 
intermediate surface is when the length scale is L = 100 
or equivalcntly T = A8K. Here, the surface is clearly 
more peaked than the original surface and the pumpi- 
ong contour encloses a reasonable fraction of the non- 
zero transmission. The pumped charge turn out to be 
0.29. The bottom surface is much more highly peaked 
and is the transmission at the length scale L = 1000 or 
at Tp = 18mK. Here, clearly, almost all the non-zero 
transmission is enclosed by the pumping contour and the 
pumped charge is 0.62. This narrowing of the transmis- 
sion peaks is expected from Luttinger liquid studies of the 
double barrier [43] , and this narrowing is what is respon- 
sible for the fact that the dissipative term in Eq. 4 starts 
contributing less and less. As shown in Fig. 3, this leads 
to charge quantisation at very long length scales or very 
low temperatures. Experimentally, it should be possible 
to study the pumped charge as a function of the tempera- 
ture and extract the interaction parameter a by plotting 
Q versus T. In Fig. 3, we have also plotted the change in 
the pumped charge as a function of the length scale for 
three different values of the interaction parameter a. As 
expected, it reaches close to quantisation earliest for the 
highest value of a. We have not taken very large values 
of a, since this approach to interactions is perturbative 
and will not work for strong inter-electron interactions. 
For strong inter-electron interactions, one needs boson- 
isation [35]. The pumped charge is perfectly quantised 
whenever the backscattering potential leads to insulating 
behaviour. The charge quantisation in a double barrier 
open quantum dot has also been attributed to Coulomb 
interactions in Rcf. [32]. The advantage of our method 
of introducing interactions perturbatively is that we can 
study the crossover from non-interacting electrons with 
low values of pumped charge to interacting electrons with 
quantised pumped charge at low temperatures. 

We can also explicitly compute how the dissipation 
term in Eq. 4 scales as we change the length scale (equiv- 
alently temperature) . The RG equation for the transmis- 
sion in Eq. 10 can be integrated to yield 

p—al f 

HI) = J ° = (19) 

where ,fo,r n are the values of the transmission and re- 
flection at I = 1, or L = d, (i.e., in the high temperature 
limit). Hence at length scales L where the second term in 
the denominator of Eq. 19 can be ignored, we find that 
t(L) cx e~ al = (L/d)~ a . Using this, we can obtain the 
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FIG. 5. Contour plots of the transmission as a function of V\ and V2 for four different lengths L = 1, 10, 100, 10,000 for the 
triple barrier system. The maxima are shown in black and the separation between the contours is 0.05, whereas the white 
regions correspond to transmission less than 0.05. The colour coding for the transmissions starting from F = 1 in black to 
F — 0.05 in light gray, is shown in the last panel. The two pumping contours are for Vb = 4.5, V p = 4.0 and 8 = tt/2 for the 
dashed contour and 8 = n/4 for the solid contour. The pumped charges for each of the contours in each of the panels is given 
in the text. The other parameters are as given in Fig. (2). 



scaling behaviour of the pumped charge as a function of 
the length scale L or the temperature T. In terms of the 
Landauer-Buttiker conductance Go = (2e 2 /^)/o, using 
Eq. 19, we obtain 

Q = Qmt --!-( — f I dtl(t) 

Aire v F Jo 

*"'^MWW- ,20) 

5 = 9 — <j) and as earlier, we have used the units Ti = ks = 
1- Qint is the integer contribution of the first term in Eq. 
4. Thus, as a function of temperature or the length scale 
L, the pumped charge scales as 

Q - Qint - c(^)- 2a = c'T 2a , (21) 
a 

at low temperatures (or long length scales), where the 
term that appears in the denominator of the integrand 
in Eq. 20 can be ignored. This can be seen in Fig. 4 
where we note that when Q — Q m t is plotted against 
L, at large values of L, the different curves fall on top of 
each other. Note that, as we mentioned before, for length 
scales above the system size L s , the scaling is no longer 
in terms of temperature, but instead is replaced by 



q - Qmt = c (^r 2a = c " L s 2a ■ (22) 

The same scaling graph is applicable here, since we have 
given it in terms of a generic L. 

• Double dot case or rib = 3 : 

Here, we have computed the 5*-matrix by first obtain- 
ing the M-matrix for 3 barriers and obtained the trans- 
mission and reflection coefficients, their phases and the 
pumped charge. The contour plots for the transmission 
as a function of the barrier heights V± and V2 is given in 
Fig. 5. The four panels which are for four different val- 
ues of the length scale, starting with the unrenormalised 
values of the transmission for L = 1, show how the peak- 
ing of the transmission occurs as we go to longer length 
scales or lower temperatures. As explained in the figure, 
the maxima (T very close to unity) are shown in black. 
The separation of transmission between the contours is 
0.05 and the lowest value of transmission (T = 0.05) 
is for the lightest of the gray contours. The new fea- 
ture that appears here is the presence of more than one 
transmission maxima for a given Ep. Note that the true 
maxima appear close to the V\ = V2 line. This is also 
true for the double barrier case (see Fig. 2). The reason 
is simply because the transmission maximum (=1) can 
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FIG. 6. Q (in units of e) as a function of the RG length scale, L for the two different contours shown in Fig. 4 (triple barrier 
system). The appropriate reflection phases for the two contours are also shown. We have set Vq = 4.5, V p = 4.0, and S — n/4 
for Ci and S = n/2 for C2. We have chosen Ef = 4.6. The other parameters are the same as in Fig. (2). 



only be reached for symmetric barriers. Here, the barri- 
ers are time- dependent and are not always symmetric, 
but the maxima still appear when the barriers are sym- 
metric. Note also that in Fig. 5, there are maxima when 
one of the barrier strengths goes to zero (or one of the 
barriers is switched off). This is simply the reflection of 
the resonance maximum of the double barrier case and is 
not relevant for the genuine triple barrier problem. The 
pumped charge is now seen to depend on the pumping 
contour chosen (i.e., the values of the barrier parameters 
that are modulated). Two possible pumping contours are 
shown in the Fig. 5, with the solid contour C\ enclos- 
ing both the transmission maxima and the dashed con- 
tour Ci only enclosing one. At high temperatures, (first 
panel in Fig. 5), the regions of non-zero transmission are 
spread out and the pumped charge can be anything since 
the dissipation term is large. It cannot be predicted and 
there is not much difference in the pumped charge of the 
two contours, which are found to be Qc\ — Qc 2 = 0.38 
for both the curves. But as we decrease the tempera- 
ture, we see that the regions of high transmission are 
squeezed together and by the time L = 100, the two max- 
ima are quite distinct. Both for L — 10 and for L = 100, 
the solid contour C\ which encloses more of the areas of 
non-zero transmission has a higher value of the pumped 
charge. For the second panel L = 10, Qc x = 0.73 and 
Qc 2 = 0.57 and for the third panel, L = 100, Qc x = 1.12 
and Qc 2 = 0.72. By the time L = 10, 000 (at low tem- 
peratures), the last panel in Fig. 5 clearly shows that the 
dashed contour C2 encloses all of the transmission around 
one of the maxima and no part of the transmission of the 
second maxima. Hence, for this contour, the the charge 
is quantised to be almost one (Qc 2 — 0.94). Ci, on the 
other hand, encloses both maxima and most of the non- 
zero transmission of the second maxima as well. Here, 



the pumped charge Qc x = 1.7, and further renormalisa- 
tion is needed for it to reach the quantisation value of 2. 
Thus, for the two barrier case, the pumped charge can 
be quantised to be one or two depending on the pump- 
ing contour chosen. This is also shown in Fig. (6), where 
we also show the change in the reflection phase for both 
contours. Clearly, the phase change is directly related to 
the first term in Eq. 4 , and hence to the charge pumped. 




FIG. 7. Contour plots of the transmission as a function of 
Vi and V2 for four barriers. The conventions are the same 
as in Fig. 5. The two pumping contours have 8 = 7r/8, 
with the solid contour Ci having Vo = 5.5 and Vp = 5.0 and 
the dashed contour C2 having Vo = 4.0 and V v = 3.5. The 
reflection phases for these contours are shown in Fig. (9). 

• Triple dot case or nt, = 4 : 

Here again, we compute the transmission, the phases 
and the pumped charge and the contour plots of the 
transmission are shown in Figs. 7 and 8. The main new 
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FIG. 8. Contour plots of the transmission as a function 
of Vi and V2 for four barrier with two pumping contours 
C3 enclosing the transmission maximum and C4 not enclos- 
ing the transmission maximum. C3 is defined by Vo = 5.5 
and V p = 2.0 and C4 is defined by V01 = 3.5, V02 = 7.0 and 
V p = 2.0. 8 — ir/2 for both contours. 

feature that appears here is the unusual shape of the res- 
onance maxima - i.e., in one case, it is almost fiat in one 
direction and appears as a line. Also, here, we find that 
there are three (non-trivial, genuine, quadruple barrier) 
transmission maxima, which appear along the V\ = V% 
line. (The reason for this has been mentioned earlier. 
Also, as has been earlier seen, the other maxima that 
occur for either V\ = or V% = and are the maxima 
through triple or double barriers.) The charge pumped 
is a function of the contour chosen. We can choose a con- 
tour that encloses all three maxima, (Ci in Fig. 7) or one 
that only encloses two transmission maxima (C2 in Fig. 
7). In Fig. (9) , the corresponding phase change in the 
reflection phase is plotted and we see that depending on 
whether the contour encloses three or two maxima, the 
phase change (and consequently, the charge pumped at 
low temperatures) is three or two units. In Fig. (8), we 
illustrate yet another important feature. The pumping 
contour has to actually enclose the centre of the reso- 
nance line; otherwise, the pumped charge is zero at low 
temoeratures. This is seen in Fig. 9, where the phase 
change for the contour that encloses the central point 
(C3) is shown to be 2tt and the phase change for the 
contour that does not (C4) is shown to be zero. 

These results can easily be extrapolated for the case 
with arbitrary number of barriers n&. The results are 
very similar. For n& barriers, there are n& — 1 true max- 
ima close to the V\ ~ V2 line. If we can choose a pumping 
contour that encloses all these maxima, the drop in the 
reflection phase will be (jib — l)27r and the maximum 
pumped charge will be (nj — l)e. To obtain quantised 
values of the charge, we need to go to low temperatures, 
so that the dissipative term vanishes. However, to obtain 
the value of the pumped charge at low temperatures, no 
explicit renormalisation or computation of transmissions 



at low temperatures need be performed. We only need to 
compute the change in the reflection phase in one pump- 
ing cycle for any given contour. This value is quantised 
and is a measure of the quantised charge that can be 
pumped at low temperatures. 

Note that all the qualitative features that we have men- 
tioned above for the pumped charge for multiple barriers 
do not depend on the details of the potentials used, such 
as the ratio of V p to Vo, or the value of Ef or the value 
of a. We choose Ep > Vi to be in the resonant tun- 
neling limit. Then as we tune the Ep for fixed values 
of the potential, the transmission amplitudes show wide 
maxima with multiple (nt, — 1) peaks, periodically as a 
function of Ep . Similar features are seen for different val- 
ues of Ep. Typically, we have chosen Ep values where 
the pumped charge is maximised. Similarly, as a function 
of the barrier separation a, the pumped charge shows pe- 
riodic behaviour, and we have chosen an arbitrary value 
of the separation. 

Note also that although these results have been ob- 
tained for quantum wells, they are also applicable to 'sin- 
gle level' quantum dots, where there is only one energy 
level (close to the Fermi energy) playing a role in the 
pumping. In other words, as long as the pumping fre- 
quency is low, so that huj < A, where A is the energy 
separation between the single level and the other levels 
in the quantum dot, the above analysis should hold. 

To conclude, in this paper, we have examined the ef- 
fects of electron-electron interactions on adiabatic quan- 
tum pumps and have shown that the limit of 'optimal' 
quantum pumps ( pumps with no dissipation) are realised 
at low temperatures, when the effects of interaction drive 
the dissipative term to zero. We have obtained the scal- 
ing of the dissipative term as a function of temperature 
and shown that it scales as T 2a vanishing at T = 0. Fu- 
ture extensions include the study of interactions on a spin 
pump [44] 
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